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ABSTRACT 

The motion of satellite galaxies around normal galaxies at distances 50-500 kpc provides a 
sensitive test for the theories. We study the surface density and the velocities of satellites around 
isolated red galaxies in the Sloan Digital Sky Survey. We find that the surface number-density of 
satellites declines with the projected distance as a power law with the slope —1.5 — 2. The rms 
velocities gradually decline: observations exclude constant velocities at a ~ 10a level. We show 
that observational data strongly favor the standard model: all three major statistics of satellites 
- the number-density profile, the line-of-sight velocity dispersion, and the distribution function 
of the velocities - agree remarkably well with the predictions of the standard cosmological model. 
Thus, that the success of the standard model extends to scales (50-500) kpc, much lower than 
what was previously considered. MOND fails on these scales for models which assume any single 
power-law number-density profile of satellites and any constant velocity anisotropy by predicting 
nearly constant rms velocities of satellites. Satellite data can be fit by fine-tuned models, which 
require (1) specific non-power-law density profile, (2) very radial orbits at large distances (velocity 
anisotropy (3 = 0.6 — 0.7 at R = 200 — 300 kpc), and (3) 2-2.5 times more stellar mass than what is 
found in the galaxies. The external gravity force - a necessary component for MOND - makes the 
situation even worse. We argue that a combination of satellite data and observational constraints 
on stellar masses make these models very problematic. 
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Introduction 



cosmology: theory — dark matter — galaxies: halos — galaxies: structure — methods: 



One hundred years after Einstein, the theory 
of general relativity (GR) is still our best theory 
of gravity. In the framework of GR, the stan- 
dard model of cosmology (AC DM) provides a 
successful description of the Universe. In this 
model, the same fluctuations which give rise to 
the observed small variations in the tempera- 
ture of the cosmic microwave background (CMB) 
grow under the force of gravity, and eventually 
form observed galaxies and other nonlinear struc- 
tures such as filaments, voids, groups and clus- 
ters of galaxies. According to the model, only 
~ 4% of the density in the Uni verse is pro- 
vided by normal baryonic matter ( Spergel et al.l 



20061 ). The AC DM model requires two ad- 
ditional components: a non-baryonic cold dark 
matter (CDM), which contributes about 30% of 
the average density of the Universe, and an even 
more mysterious dark en ergy, which makes up the 
rest ( Spergel et alJlioOrt ). The model is remark- 
ably successful on scales larger than a few Mega- 
parsecs. It predicted the amplitude and the spec- 
trum of angular fluctuation s in the CMB and in 



the distribution of galaxies (jBardeen et al. 111987 



Holtzman 19891) that were later confirmed by ob- 



serva t ions ( Spergel et al.l 120061; Netterfield et al. 



20021: ICole et ah I 120051: iTegmark et al. I 12004) 



However, the ACDM model faces challenges on 
smaller scales. The most difficult ones are re- 
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lated with the rotation in the inner parts of spiral 
galaxies. It seems that the theory predicts too 
much dark matter inside ~ lkpc from the cen- 



ters of galaxies (IMoore Ill994t iFlores fc Primack 



1994 Ide Blok et al. 112001ft . While th ere are some 
possible solutions of the problem dRhee et al. 



2004; H avashi fc Navarro H2006tlValenzuela et al. 
2007), the problems on small scales are the 
strongest challenge the standard model has en- 
countered. When compounded with the fact that 
there is no direct evidence of dark matter or dark 
energy, the current problems of the standard cos- 
mological model have encouraged a small but 
growing community of physicists to propose al- 
ternative theories of gravity to avoid the need for 
dark matter. 

This is the case for Modif ied Newtonian D y- 
namics (MOND), proposed bv lMilgronil (|l983ft to 
explain the rotation of galaxies without dark mat- 
ter. According to MOND, the rotation curves in 
the outer regions of galaxies do not decline because 
the force of gravity is significantly stronger than 
for Newtonian gravity. At early times MOND's 
main appeal was its simplicity: there is no need 
to make the assumption that the Universe is filled 
with particles that nobody has seen. Additional 
motivation came later from difficulties with ex- 
plaining anomalies in the trajec tories of the Pi- 



oneer 10 and 11 space missions ([Andreson et al. 
2002ft . 



Yet, for a long time MOND was not more than a 
conjecture. Only recently, Bekenstein proposed a 
relativistic ve rsion named tenso r vector scalar the- 
ory (TeVeS) (|Bekenstein 112004ft . This alternative 
theory of gravity provides a framework to make 
predictions of numerous important observational 
phenomena, which ACDM has already done: the 
temperature fluctuations seen in the CMB, grav- 
itational lensing, and the large scale structure of 
the universe. With maturity came problems. Ro- 
tation curves of some galaxies - the initial strong 
argument for MOND - cannot be explained by 
MOND. In about 1/4 of galaxies considered by 
proponents of MOND the predicted velocities are 
well abo ve the observations in the very central 



regions ([Sanders fc McGaugh 1 12002ft . RMS ve- 



locit ies of stars in some dwarf sph eroidal galax- 
ies (|Lokas. Mamon. fc Prada 1 12006ft also present 
problems. 

So far, the most severe challenges for MOND 



are coming from clusters of galaxies. Dynam- 
ics of galaxies in clusters cannot be explained by 
MOND and requires introduction of dark matter, 
possib ly in the form of a massive (~ 2eV) neu- 



trino ([Sanders fc McGaugh I [2002ft . We do not 



know whether this modification can explain some 
properties of clusters of galaxies such as the "Bul- 
let Cluster", where the baryonic mass (galaxies 
and gas) is clearly separated from the gravita- 
tion al mass, as indicated by gravitational lens- 
ing ( Clowe et al. 1 120061: lAngus. Famaev. fc Zhao 



2006). In any case, for MOND to survive too it 



must invoke dark matter in the form of massive 
neutrinos and dark energy in the form of an ar- 
bitrary constant added to a combin ation of two 



scala r fields used in TeVeS MOND ([Bekenstein 
l2004ft . 

There is no doubt that alternative theories 
of gravity represent a challenge to the standard 
model of cosmology and GR. Any theory or model 
must respond to these challenges. Here we present 
a number of observations to test gravity and dark 
matter in the peripheral parts of galaxies at dis- 
tances 50-500 kpc from the centers of galaxies. 
These scales can be tested by studying the mo- 
tion of satellites of galaxies. This is a rela- 
tively old field in extragalactic astronomy and his- 
torically it was one of main arguments for the 
presence of dark m atter ([Zaritskv fc White lfl99l 
IPrada et aL~ll2003ft . 

The paper is organized as follows. In Section 2, 
we present the observational results drawn from 
the SDSS and the predictions from the standard 
model of cosmology. Predictions from MOND are 
computed and discussed in Section 3. Finally, con- 
clusions are given in Section 4. 



2. Observational results 

We use the Sloan Digital Sky Survey (SDSS; 
www.sdss.org) - the largest photometric and 
spectroscopic astronomical survey ever under- 
taken of the local Universe - to study the mo- 
tion of satellites. As of Data Rele ase Four (DR4) 
( Adelman-McCarthv et al. 2006ft . imaging data 
are available over 6670 deg 2 in five photomet- 
ric bands. In addition to the CCD imaging, the 
SDSS 2.5m telescope on Apache Point, New Mex- 
ico, measured spectra of galaxies, providing dis- 
tance determinations. Approximately half mil- 
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Fig. 1 . — Left panel: RMS of the linc-of-sight velocities of satellites orbiting red isolated galaxies with absolute 
magnitudes indicated in the plot. The luminosity decreases from top to the bottom curve. Satellites move 
faster around more luminous galaxies. In all samples of galaxies the rms velocity is declining with distance. 
Right top panel: Distribution of observed line-of-sight velocities (dots with error bars) has simple structure of 
the Gaussian distribution with a small constant background (full curve). The background is due to objects, 
which in projection lie close to the central galaxy, but that are far from it in 3D space. The dashed curve is 
the distribution of velocities expected in the standard cosmological model drawn from AC DM simulations. 
Right bottom panel: Surface number-density of satellites (dots with error-bars; arbitrary units) orbiting 
galaxies with luminosity —20.0 < M g < —21.5. The dashed line is the power law with a slope —2. The full 
curve is the prediction of the ACDM model. 



lion of galaxies brighter than r = 17 '.77 over 
4700 deg 2 have been targeted for spectroscopic 



observations as part of SDSS and are included in 
DR4. Rcdshift accuracy is better than 30 km/s 
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and the overall completeness is ~90%. For our 
study we compute rest frame absolute magni- 
tudes in the g-band from the extinction-corrected 
apparent magnitudes assuming a ACDM cos- 
mology with a Hubble constant h = 0.7 (Hq = 
100/ikms" 1 Mpc _1 ). Galaxies are split into red 
(early-types) and blue (late-types) populations 
based on the bim odality obseryed in the u — r 



color distribution (Ba ldrv et al. I l2004l . e.g.). The 



local minima between the peaks of the color distri- 
bution occur near u — r = 2.3. All magnitudes and 
colors are k-corrected to z = 0. Because calcula- 
tions of MOND gravity for non-spherical objects 
are complicated, we restrict our analysis only to 
red galaxies, the vast majority of which are either 
elliptical galaxies or are dominated by bulges. Our 
galaxy sample was selected from the full redshift 
sample by taking all galaxies with recession ve- 
locity 3000 km/s < cz < 25000 km/s. The to- 
tal number of selected galaxies is about 215,000. 
The SDSS heliocentric velocities were converted 
to the Local Group standard of rest before com- 
puting distances. We select our host galaxies as 
galaxies with absolute g-band magnitude brighter 
than M g = —19.0 and isolated: a galaxy must be 
at least 4 times brighter than any other galaxy 
within a projected distance R < 1 Mpc and a 
line-of-sight velocity difference AV < 1500 km/s. 
We define satellites as all galaxies being at least 
4 times fainter than their hosts and found within 
a projected distance R < 1000 kpc and velocity 
difference AV < 1500 km/s with respect to their 
hosts. Typically we find about 1.5 satellites per 
host. In total we have 9500 satellites with a mean 



luminosity of about M„ 



48.0. We bin the 



host galaxies by luminosity and collect informa- 
tion about the distribution of relative velocities 
AV and the number of satellites as the function 
of projected distance R. 

Figure 1 presents the observational results for 
primaries, which are at least 6 times brighter than 
any satellite. We also used primaries, which are 4 
times brighter, and primaries, which arc 10 times 
brighter than their satellites. We do not find any 
trend with the primary-satellite magnitude gap: 
results are nearly the same (see Section 3 for de- 
tails). The only difference is the statistics, and, 
thus, the error-bars. The distribution of line-of- 
sight velocities (top right panel) clearly shows a 
two-component structure: a homogeneous back- 



ground of interloper galaxies (dwarfs which hap- 
pen to lie along the line-of-sight, with large physi- 
cal distances from the host galaxies but small pro- 
jected and velocity differences, but which are not 
associated with the host) and a nearly Gaussian 
component. The surface density of the satellites 
also shows the same structure: at large separa- 
tions the number density goes to a constant due 
to interlopers. We subtract the background and 
plot the surface density of satellites in the bottom 
right panel. In order to study the velocity dis- 
persion of satellites AV^ = (AV 2 ) as a function 
of projected distance R we use a maximum likeli- 
hood method to approximate the number of satel- 
lites N (R, AV) per unit projected area and per 
unit velocity difference AV using an 9-parameter 
function in the form of a constant plus a Gaussian 
distribution with variable velocity dispersion and 
normalization: 

N(R, AV) = n(i?)exp [-V 2 /2a 2 (R)] /V2^cr+n Q , 

(I) 

where the surface number-density of satellites 
n(R) and the rms line-of-sight velocity a = 
AV lms {R) are 3th order shifted Chebyshev polyno- 
mials of the first kind. The parameter uq is a con- 
stant representing the background of interlopers. 
The left panel in Figure 1 shows the resulting rms 
velocity of satellites for three magnitude bins. In 
order to estimate statistical uncertainties, we run 
Monte Carlo simulations using the same number 
of hosts and satellites as in the corresponding mag- 
nitude bin in Figure 1. Note that the data points 
are correlated. We use Monte Carlo simulations 
to test statistical significance that the observed 
velocities are declining with distance: data for 
each magnitude bin reject a constant AV rms (i?) at 
about 3cr confidence level. We also studied blue 
galaxies and find declining velocities for them. 
These re sults are in agreement with previous es - 



timates (jPrada et al. 1120031 ; IConrov et al. 1 12007). 



but we now can exclude constant rms velocities 
much more reliably. 

3. Predictions of the ACD A/model 

In order to compare observational results with 
the ACDM predictions, we use h igh-resolution 



cosmological N-body simulations (jPrada et al 
2006), which we then treat as if they are the 
observational data. The simulations were done 
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Fig. 2. — The line-of-sight velocities of galaxies in two luminosity ranges: —20.5 < M g < —21.1 (2300 
satellites) and —21.1 < M g < —21.6 (2700 satellites). The vertical lines indicate 68% confidence levels. 
The left panel shows predictions (full curves) of the standard cosmological model for galaxies hosted by 
dark matter halos with maximum circular velocities of ~ 340 km/s (top curve) and <~ 270 km/s (bottom 
curve). The right panel shows MOND predictions. The two bottom curves are for stellar masses 8 x 1O 1O M 
estimated for galaxies with M g = [—20.5,-21.1]. The full curve is for models with isotropic velocities 
= and with the observed slope a = — 3 of the 3D number density of the satellites. The curve was so 
much below the observed data points that for other models we decided to use the slope a = —2.5, which is 
marginally compatible with the observations. The dashed curve is for orbits, which are preferentially radial 
(0 = const = 0.5). The two top curves are for the stellar mass 2 x 10 n M Q and for = (full curve) and 
= 0.6 (dashed curve). The KCDM makes quite reasonable predictions, while MOND has problems. 
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for the standard AC DAI model with parameters 
sigma$ = 0.9, flo = 0.3, h = 0.7, n = 1.. 
We use two simulations. Both simulation had 
512 3 particles. One simulation had the com- 
putational box 120 h Mpc box, mass resolu- 
tion mi — 1-07 x 10 9 /i _1 Mq and force resolu- 
tion 1.8 h~ x kpc. The second simulation has the 
computational box 80 /i _1 Mpc box, mass reso- 
lution mi = 3 x 10 h Mq and force resolution 
1.2 h~ 1 kpc. We do not find any differences be- 
tween the simulations and decided to present re- 
sults of the larger simulation, which provides bet- 
ter statistics of halos. We select isolated halos as 
halos, which within projected distance of 715 kpc 
and within relative velocity difference 1000 km/s 
do not have halos or subhalos with maximum cir- 
cular larger than 1 /2 of the halo's circular velocity. 
This corresponds to the mass ratio of a factor ten 
for distrinct halos. Circular velocities are better 
characteristics of subhalos, which have poorly de- 
fined masses. Although we resolve subhalos, we 
dicided to use dark matter particles as proxies for 
satellites. At distances larger than a fraction (1/4- 
1/3) of a virial radius from center of a halo satel- 
lites trace the motion a nd the spatial dis tribution 
of the dark matter (e.g. ISales et al.ll2007l ). 

Results presented in Figures 1 and 2 show that 
all three characteristics of observed satellites are 
reproduced by the model. Note that we actually 
do not make fits. Results from simulations have 
only one parameter: the maximum circular veloc- 
ity of the dark matter halo. We just plot what we 
obtain from the simulations. 

Once we fix the maximum circular velocity, the 
results are fairly insensitive to parameters of the 
cosmological model. Two factors may affect the 
shape of the theoretical V-R diagram: the halo 
concentration and the velocity anysotropy. The 
halo concentrations c = R v i v /R s are in the range 
c = 8 - 10 for halos (M vlr = (2 - 10) x 1O 13 M ) 
and normalization <7s cosidered in this paper 
([Bullock et alj|200ll iNeto et al.ll2007i) . (Here R vir 
is the virial radius and R s is the characteristic 
radius of the NFW profile). For virial radii in 
the range 300 — 500 kpc, the exact value of the 
halo concentration affects only the central region 
R < 50 — 70 kpc. The average halo concentration 
depends weakly on the amplitude of the pertur- 
bations erg . Existing observational data leave only 
very narrow range for variation of as = 0.8 — 0.9, 



which also limits the range of concentrations. 
The velocity anisotropy is less certain, but it is 
not large. N-body simulations indicate that (3 is 
slightly positive and has a tendency to increase 
with the distance for radii smaller than the virial 



radius (e.g. lCuesta et al.ll20071 ) with typical values 
j3 s=s 0.2 — 0.3 for halo masses considered in this 
pape r. At larger distanc es (3 declines and goes to 



zero ( Cuesta et all 2007) 



In order to demonstrate that our results of N- 
body simulations are robust and reliable, we use 
solutions of the Jeans equation, which we obtain 
for parameters compatible with numerous previ- 
ous simulations. Once a solution for the radial ve- 
locity dispersion is obtained, we integrate it along 
a line of sight with appropriate corrections for 
the velocity anisotropy. For the density use the 
NFW profiles with concentration c = 9 — 10. For 
the sake of completeness, we add the cosmolog- 
ical background density, which has only a small 
effect at large distances. For j3 we use approxi- 
mation given by eq. ([7]) with parameters chosen 
in such a way that /3(0) = and at R — 300 kpc 
j3 = 0.2—0.4. We use the same range of circular ve- 
locities as in Figure O For more massive halos we 
use Af vir = 8.5 x 10 12 h^ 1 M Q and c = 9. For less 
massive halos we adopted M v - lr = 3.8 x 10 12 /i _1 Mq 
and c = 10. Results presented in Figure [3] clearly 
indicate that very simple stationary models can 
accurately reproduce results of simulations. 

4. Predictions from MOND 
4.1. Analytics 

The situation is different for MOND because 
there are no predictions on the same level of so- 
phistication as for ACDM. In principle, those 
predictions can be made, but at this moment they 
have not been yet made. Thus, we have only one 
option: solve the Jeans equation for spherical sys- 
tems. When doing so, we have a freedom of choos- 
ing two functions: the number-density profile of 
satellites and the velocity anisotropy (3(r). We 
also have two free parameters: stellar mass M* 
and the magnitude of external force (see below for 
details). The predictions are constrained by two 
functions n{R) and AV(R). The velocity distri- 
bution function also gives constraints, but those 
are relatively weak. Only the models with a large 
\B\ > 0.85 can be excluded. The stellar luminos- 
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Fig. 3. — Comparison of results of N-body sim- 
ulations (full curves) with solutions of the Jeans 
equation (dashed and dot-dashed curves. In the 
top part of the plot the vertical error bars show 
observational constraints for galaxies in the lumi- 
nosity range —21.1 < M g < —21.6. The corre- 
sponding full curve is for halos in simulations with 
maximum circular velocity ~ 340 km/s. The two 
other curves are solutions of stationary spherically 
symmetric Jeans equation with the same maxi- 
mum circular velocity. The dashed curve is for 
slightly more radial velocities (velocity anisotropy 
(3 = 0.34 at R = 300 kpc. The dot-dashed curve 
of for a model with = 0.17 at the same distance. 
The bottom curves are for halos with maximum 
circular velocity ~ 270 km/s. The vertical er- 
ror bars are for galaxies in the luminosity range 
-20.5 < M g < -21.1. The other curves are la- 
beled in the same way as the top curves. 



ity and colors constrain the stellar mass. Roughly 
speaking we have two arbitrary functions to fit two 
observed functions. It should not be difficult given 
that the model is viable. 

In our case the Jeans equation for the radial ve- 
locity dispersion a 2 (r) can be written in the form: 



dr 



. (2/3 + a) 



-g{f),ot = \np/dhxr, (2) 



where (3(r) = l — a 2 _/2a 2 is the velocity anisotropy 
and g(r) is the gravitational acceleration. The for- 
mal solution of the Jeans equation can be written 
in the form: 



1 



X(r)p(r 
X(r) = cxp 



X(r)p(r)g(r)dr, (3) 
(4) 



It is convenient to chose forms of a and /3 such 
that the integral in eq.Q is taken in elementary 
functions. We used the following approximations: 



/3 = /3o + 



p = r- a (l + —)- a \ 
r n 



/3 = /3 , X = r Wo , 
Pi 



1 + ^' 



X = r 2/3a 



-| 2/3! 



(5) 
(6) 

1 

(7) 



/3 = A)+ 



l + (^) 2 '' 



„2/3 



r 


2j8i 


1 






l + (^) 2 



(8) 

In the case of Newtonian gravity the acceler- 
ation g = gN{r) = GM(r)/r 2 , where the mass 
M(r) includes both normal baryonic mass and 
dark matter. For MOND the acceleration of spher- 
ical systems (^mond is given by the solution of the 
non-linear equation 

5MONDM(|g M OND + gext|/ao) = GM*(r)/r 2 , (9) 



where M* is the mass of only baryons and an = 
1.2 x 10~ 8 cm sec -2 . The term g cxt is the ex- 
ternal constant gravitational acceleration. It for- 
mally breaks the spherical symmetry. Thus, the 
eq. © is only an approximation valid for small 
external force. Her e we use the same app r oxima - 
tion as eq. (10) in ISanders k. McGaugh I (|2002h . 
which we average over angles between a constant 
external acceleration g oxt and internal radial ac- 
celeration (7MOND- 
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Fig. 4. — The line-of-sight velocities of galaxies with —20.1 < M g < —20.9 (left panel) and —21.1 < 
M g < —21.6 (right panel) are compared with MOND models with carefully tuned parameters. We use 
M* = 8 x 1O 1O M0 and M* = 3 x 1O 11 M0 for the two bins. No external gravity force was assumed. Models of 
MOND can fit the data only when very radial orbits with /3o = 1.0 (eq.(7)) are used (dashed curves). Models 
with more reasonable /3o = 0.5 fail (full curves). Data for more luminous galaxies (right panel) require twice 
as much of stellar mass as actually observed in the galaxies. 



The function (J,(x) can have differe nt shapes. 



We t ried the originally proposed form (jMilgrom 
19831 ) n(x) — x/y/l + x 2 , but accepted and used 



for all our analysis the function fi(x) = x/(l + x), 
which gives slightly better results. The function 
a(r) is limited by the observations presented in 
Figure 1. The velocity anisotropy /3(r) is a free 
function, but there are constraints. Asymptoti- 
cally j3 goes to zero at large distances where grav- 
itational effect of the central galaxy diminishes. 
Ideally the velocity anisotropy at small distances 
also should be declining. There are different ar- 
guments why this should be the case. (1) The 
tangential velocities of few satellites of the Milky 
Way, for which the proper moti ons are measured, 
strongly reject radial orbits ( Kallivavalil et al.l 
120061 (Piatek et al.l [20071) . There is no reason to 
believe that our Galaxy should be special in this 
respect. (2) Experience with gravitational dynam- 
ical systems indicate that in dynamically relaxed 
systems « 0. Numerous simulations of cosmo- 
logical models illustrate this: (3 is small in the 
central regio n and increases to 0.2 — 0.3 at the 
viral radius ( Cuesta et al.l 120071 ). Note that we 



should distinguish the velocity anisotropy and the 
orbital eccentricity. For centrally concentrated 
objects, which we are dealing with here, already 
isotropic velocities imply typical peri-/apocenter 
ratios of 1:4-1:5. If eccentricity is larger, a signif- 
icant fraction of satellites comes too close to the 
central ~ 10 kpc region where the satellites are 
destroyed by tidal forces. 

The term g ext in eq.© describes the effect 
of external gravitational field. It is specific for 
MOND. In Newtonian dynamics a homogeneous 
external gravity does not affect relative motions 
inside the object. Because the MOND gravity is 
nonlinear, the internal force is affected by the ex- 
ternal field (note that this is not the tidal force). 
This external effect is quite complicated. The 
magnitude of the external f orce is substantial fo r 
the motion of th e satellites. iFamaey et alJ (|2007l ): 



Wu et all (|2007l ) point out that in MOND numer- 



ous sources and effects generate about equal mag- 
nitude of gext- For example, 600 km/s motion of 
the Local Group relative to the CMB implies ac- 
celeration of 600 km/s/10 Gyrs= 0.015ao- Here we 
assume that the acceleration is constant over the 
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whole age of the Universe. If g ext increases with 
time, as it may be expected, then the acceleration 
is even larger. Infall with 250 km/s in the direc- 
tion of Virgo gives about half of the value. M31 
produces about the same magnitude of the accel- 
eration. While the acceleration in MOND does 
not add linearly, it is reasonable to assume that 
9rmext > 0.015ao- We will explore the effect the 
external field later. 

4.2. Comparison of MOND models with 
observations 

In order to make MOND predictions, we must 
estimate the stellar mass for galaxies in our anal- 
ysis. For two subsamples presented in Figure [4] we 
use the magnitude bins M g = -20.5 - 21.1 (2400 
satellites) and M g = -21.1-21.6 (2700 satellites). 
The average luminosities of galaxies in the bins are 
L g = 2.0 x 10 10 i Q and L g = 4.0 x 10 10 L Q . When 
estim ating the luminosit ies, we assume M S!0 = 
5.07 (jBlanton et al.ll200lh . Using measured u — g 
colo rs we estimate the s tellar masses of galax- 



ies dBell fc de .Tonal l200lh : M* = 7.2 x 1O 1O M 



and M* = 1.5 x 10 11 M Q for the two luminos- 
ity bins. This implies the mass-to-light ratios are 
nearly the same for the bins: M/L g = 3.7. These 
M/L estimates are close to predictions of stel- 
lar population models. In our case, the galaxies 
are red (u — g > 2.3), old mostly ellipticals, for 
which we expect nearly solar metallicity. Indeed, 
Marastonl ( 2005 ) gives M/Lb = 5 for the Kroupa 
IMF for stellar population 10 Gyrs old. Adjust- 
ing for the 0.4 mag difference between B and g 
bands, we get M/L g = 3.5 (the Salpeter IMF gives 
M/L g — 5.5). Larger metallicity is highly unlikely 
(the galaxies are not really massive ellipticals) and 
it does not make much difference: 10 percent in- 
crease if we take twice the solar metallicity. 

Using different functional forms for j3(r) , we 
solve the equations numerically and then integrate 
the solution along the line-of-sight to get a pre- 
diction for AV rms . Figure 2 presents results for 
the two different magnitude ranges and for differ- 
ent parameters (3 =const and a =const. MOND 
definitely has problems fitting the data because 
for any constant (3 it predicts a sharply declining 
rm velocity at small distances followed by nearly 
fiat velocities at large distances. The data show 
just the opposite behavior. It is easy to under- 
stand why predictions of MOND has this shape. 




Fig. 5. — RMS velocities of satellites for galaxies 
selected by stellar mass with the average M* = 
2.2 x 10 n M Q . The top full curve with the er- 
ror bars is for a sample of 2600 satellites, which 
are more than 6 times dimmer than the primary 
galaxy. The dashed curve is for more stringent 
isolation criterion of satellites 10 times less bright 
than the primary (995 satellites). The differences 
are not statistically significant and are less than 
10 km/s for R > 50 kpc. The bottom curve is for 
MOND with the same stellar mass as in the obser- 
vations and with tuned parameters to produce the 
best fit. The MOND prediction falls much below 
the observational data. 

At large distances the Newtonian acceleration is 
very small and the acceleration is strongly dom- 
inated by the MOND correction: (?mond oc 1/r. 
In this case, the solution of eq.(l) gives a =const 
for any constant value (3 and a. 

We can try to salvage MOND, but so far th ere is 
no simple way of doing it. Angus et al. ( 20071 ) find 
that the way to improve the situation is to have 
a model with variable slope a and with extremely 
radial orbits. We also tried different functional 
forms and numerous combinations of parameters. 
While the parameters, which fit the surface den- 
sity and the rms velocities vary, they all give the 
same answer: density slope in the central 40 kpc 
should be a > —2.0 and the velocity anisotropy 
at large distances should be (3 > 0.7. Figure H] 
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Fig. 6. — The ratio of mass required by MOND 
to the stellar mass for galaxies selected by stellar 
mass. For small mass galaxies MOND requires 
masses, which are in reasonable agreement with 
the observed stellar masses. Yet, for more massive 
galaxies MOND needs more mass than actually is 
observed. 

shows (dashed curve s) the results for the same set 
of parameters as in Angus et al.1 ( 20071 ): a = 0, 
ai = 3.1, r a = 40 kpc, r b = (20 - 40) kpc, f3 = 1 
and Pi = —2.57. Here we use the approximation 
given in cq.([7|). The solution is very contrived: 
relatively small deviations from the best behavior 
(small slope in the center and radial orbits in the 
outer radii) result in failed fits. The full curves in 
the Figure show what happens when the velocity 
anisotropy gets less radial: f3 = 0.5. 

Even this fine-tuned solution fails unless the 
mass-to-light ratio for the larger magnitude bin 
is arificially increased by a factor of t wo: M/L a = 
7.5, w hich gives M* = 3 x 10 n M Q . lAngus et a! 



([20071 ) also found the same trend, but they made 
two mistakes, which did not allow them to clearly 
see the problem. First, the solar mass-to-light ra- 
tio was used for the B band instead of the g band. 
Second, the width of the magnitude bin is pre- 
sented as an error in M/L giving impression of 
very large uncertainties in M/L. This is not cor- 
rect: the statistical uncertainty of the average lu- 
minosity of galaxies in each bin is very low and 



can be neglected. It should be noted that there is 
nothing special about the galaxies, which are used 
here. The average luminosities differ by a factor 
1.7. So, it is not a large difference. Colors of 
the galaxies are practically the same, which then 
gives the same M/L if we use stellar population 
models. The rms velocities are also perfectly con- 
sistent with simple scaling. For example, the ratio 
of rms velocities at the same projected distance of 
70 kpc is 1.3 implying simple scaling L cx a 2 R. 
Roughly speaking, we double the luminosity and 
that doubles the stellar mass. This does not work 
for MOND: it needs twice more stellar mass. This 
is definitely a problem because there is no justifi- 
cation why galaxies with the same colors, with the 
same old population and practically the same lu- 
minosity should have dramatically different IMF. 
The differences are very large: most of the stellar 
mass in the more luminous bin should be locked up 
in dwarfs with ~ 0.2M©, while there is relatively 
little of those in the lower bin, which is consistent 
with the Kroupa IMF. 

In order to make the argument even more clear, 
we make analysis of velocities in a different way. 
This time we split the sample by stellar mass, but 
we still keep only red primaries with u — r > 2.3. 
We make the analysis twice: for satellites 4 times 
and for satellites 10 times less bright than the pri- 
mary galaxy. There is no systematic difference 
between the two isolation conditions: within la 
the results are the same. Figure [5] illustrates this 
point. We select primary galaxies to have the stel- 
lar mass in the range M* = (1.6 - 3.2)10 n M©. 
The average stellar mass is (M*) = 2.2 x 10 n M©, 
and the average luminosity is (L g ) = 5.4 x lO lo i0. 
There is a hint that more isolated primaries have 
slightly larger velocities of satellites. Still, the dif- 
ferences are not statistically significant: for radii 
larger than 50 kpc the differences are smaller 
than 10 km/s. Then we take the observed stellar 
mass and use it for MOND models and apply the 
best tuned p arameters. W e find that parameters 
suggested by lAngus et al. ( 2007 ) (r^ = 40 kpc, 
Po = 1, Pi = —2.57) improve the fits as compared 
with a constant P models. Still, they are not ac- 
ceptable. For example, at R = 150 kpc the MOND 
model is a 7a deviation. We found a better solu- 
tion (rj, = 30 kpc and other parameters the same) , 
which places MOND "only" at 6.2a. 

In order to envestigate what stellar mass is 
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needed for MOND, we split the sample of red 
primaries into 5 mass bins ranging from M* ~ 
5 x 10 10 A/ Q to M« ~ 5 x 1O U M . We chose 
the less stringent isolation condition because it 
gives typically 2-2.5 time more satellites result- 
ing in smaller statistical errors. Each bin has a 
large number of satellites: 1200 — 2900. We then 
run a grid of MOND models with different masses 
and select models, which make best fits. Just as 
the MOND models in Figure [H there are models, 
which marginally fit the data. For each model we 
get stellar mass required by MOND and compare 
it with the stellar mass estimated by the stellar 
population models. Results are presented in Fig- 
ure O For smaller masses MOND gives masses, 
which are compatible with actual stellar mass ob- 
served in the galaxies. The plot clearly shows the 
problem with massive galaxies: MOND requires 
increasingly more mass than observed in the galax- 
ies ending up in large (~ 2.5) disagreement with 
the observations. 

The external gravity force g cx t is another 
MOND component. It must exist on the level 
of 5oxt ~ O.Olao- Figure [7] shows what happens to 
the models when we add the external force. As the 
starting models we use the best fits shown in the 
Figure [4] as dashed curves. The models with the 
external force make the fits much worse. Taken at 
face value, the models with realistic g cx t — 0.015ao 
can be rejected. Again, we can make the model 
work if we increase M/L by a factor ~ 1.5. Yet, 
this will make the situation with stellar masses 
even worse, than what we already have. 

5. Discussion and Conclusions 

Using the SDSS DR4 data we study the distri- 
bution and velocities of satellites orbiting red iso- 
lated galaxies. We find that the surface number- 
density of the satellites declines almost as a power 
law with the slope —2.5 — 3. The distribution of 
the line-of-sight velocities is nearly a perfect Gaus- 
sian distribution with a constant component due 
to interlopers. The rms velocities are found to 
gradually decline with the projected distance. The 
constant rms velocity (isothermal solution) can be 
rejected at a 10 a level. 

Observational data strongly favor the standard 
cosmological model: all three major statistics of 
satellites - the number-density profile, the line- 



of-sight velocity dispersion, and the distribution 
function of the velocities - agree remarkably well 
with the theoretical predictions. Thus, the suc- 
cess of the standard model extends to scales (50- 
500) kpc, much lower than what was previously 
considered. 

MOND fails badly in cases with any acceptable 
power- law approximation for the number-density 
of satellites and a constant velocity anisotropy by 
producing sharply declining velocities at small dis- 
tances followed by nearly flat velocities at large 
distances - just the opposite of what is observed 
in real galaxies. Models may be made to fit the 
satellites data only when all the following condi- 
tions are fulfilled: 

1. The slope of the density in the central 50 kpc 
is less then -2. 

2. Satellite velocities are nearly radial in outer 
regions (/3 « 0.6-0.7 at R = 200-300 kpc). 

3. Mass-to-light ratio increases from M/L g w 



3.5 at M g 
17) at M g 



-20.9 to M/Lg 
-22. 
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4. External force of gravity is smaller than the 
expected value of g ox t = 0.015ao. The neg- 
ative effect of the external force can be re- 
moved by additionally increasing the M/L 
ratios. 

We find that the later three conditions are dif- 
ficult to realize in nature. Satellites do not fall in 
to their parent galaxies with zero tangential veloc- 
ities. The velocities are induced by other neigh- 
boring galaxies and by large-scale structures such 
as filaments. To some degree, it is similar to the 
tidal torque, which is responsible for the origin 
of the angular momentum of galaxies. Yet, the 
interactions are more efficient for providing ran- 
dom velocities. Measurements of peculiar veloci- 
ties of galaxies the Local Volume (3-5 Mpc around 
the Milky Way) found deviati ons from the Hubbl e 



flow that about 70-80 km/s (jMaccio et al.ll2005h . 



At these distances the gravitational pull of other 
large galaxies in the aria is larger than that of the 
Milky Way and we expect that the same deviations 
exist for the perpendicular velocity components. 
When a satellite with 70-80 km/s falls falls from 
1 Mpc and gets to 400 kpc from MW its tangen- 
tial velocity increases few times (conservation of 



11 




80 



100 



200 
R(kpc) 



300 



100 



100 



200 
R(kpc) 



300 



400 



Fig. 7. — Effects of external gravity force. Parameters are the same as in Figure 3. For both low (left panel) 
and high (right panel) luminosity bins we assume the most favorable MOND parameters: radial orbits and 
adjusted M/L (models shown with dashed curves in Figured]). The magnitude of the external field is shown 
in the left panel. Models cannot fit the data without additional increase of M/L. 



the angular momentum). So, we expect (3 < 0.5. 
We have the same measurements for the satellites 
in SDSS. When we look at distances 0.8-1 Mpc, we 
find that the rms velocities hardly correlate with 
the luminosity and are about 100 km/s. Note that 
AC-DMperfectly fits the constraint: at these dis- 
tances it predicts j3 < 0.4. Thus, MOND assump- 
tion of very radial orbits does not agree with the 
observations. 

The large M/L ratios is another problem for 
MOND. So far M/L was treated almos t as a 
free parameter ( Sanders k, Noordermeerl 12007 : 



ity. IZoccali et all (|2003f) give the M/L v = 3.67, 



Angus et al.1 120071 ). This is not correct. It 



should be reminded that observationally there 
are no indications of large variability of IMF for 
vastly different stellar ages and metallicities. The 
IMF shows similar flattening for masses smaller 
than Mq in the solar ne ighborhood and in the 



spheroid (IChabrie: 
(IZoccali et alJl200, 



200 lh . in the galactic bulge 



and in the Ur s a Min or dwarf 
spheroidal galaxy (jFeltzing et al. 1999f ). The 
galactic bulge has parameters similar to the red 
galaxies studied in this paper: very old (on aver- 
age) stellar population with nearly solar metallic- 



which corresponds to M/L g — 3. This is close to 
what we found for the red galaxies in SDSS. 

In the case studied here, MOND requires that 
M/L changes by a factor 2-2.5 when there is 
no change in colors and the absolute magnitude 
changes only by AM w 0.5. There is no justifica- 
tion for this to happen. The situation is similar 
to what ISanders fc Noordermeer ( 20071 ) found for 
early type large spirals where MOND models fail 
for any IMF discussed in recent 20 years. Only 
the Salpeter IMF gives better results. Even in this 
case the models deviate by a factor of two around 
the stellar population predictions. The combina- 
tion of all these problems makes MOND a very 
implausible solution for the observational data on 
satellites of galaxies. 

There is one significant difference between 
MOND and AC DM. The latter makes predic- 
tions for the distribution of mass and velocities 
for isolated galaxies and those predictions match 
the observational data. There are no theoreti- 
cal predictions for MOND. What we conveniently 
called "MOND predictions" were actually require- 
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merits. For example, MOND must produce the 
radial velocities of satellites or shallow slope of 
the number-density in the central region. 
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